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Abstract 

This paper deals with several ways to generalize the Flanders’ theorem to matrix triples. We 
consider six invertible matrices and try to write them as the possible products of three matrices. 
Initially, we describe a wide set of necessary conditions so that this system be solvable, show¬ 
ing that they are not sufficient. Next, we study the simultaneous solvability of two equations, 
selected appropriately among the matrix system. The rest of the paper is devoted to the study 
of a particular case, in which the six given matrices are simultaneously diagonalizable, with 
distinct nonzero eigenvalues. In this case, we obtain a necessary and sufficient condition for 
the solvability of the full matrix system. Moreover, an explicit solution to it is constructed. 
Certain technical results necessary for this work may be of independent interest. 
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In [1] it was first observed which pairs of matrices B\ e M m (C) and Bj e M n (C) 
may be written as B\ = A\Ay and By = AyA\ with A i e M mn (€) and A 2 e 
M n , m (, C). Subsequent work [8] has refined and given a straightforward proof [6] 
of this fundamental result. In the event that m — n and B\ and By be invertible, 
it is necessary and sufficient that B\ and By be similar, which, by itself, is 
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easily proven. Recently, there has been interest in similar characterization of products 
among k ^ 3 matrices. Here, because the number of possible products, using each of 
the k matrices exactly once, grows factorially in k, several analogous questions are 
possible. For example, in [ 2 , 3 ], the ^-tuples B\, .... B k (not necessarily invertible) 
which may be written as cyclic products among A/,..., A k , not necessarily square, 

(e.g., B\ = Ai A 2 ■ • • A k , B 2 = A 2 • • ■ A k A\, . .., B k = A k A\ • • ■ Aa-_i) were char¬ 

acterized. Again, the square, invertible case is relatively straightforward. Here, we 
consider two other natural generalizations and concentrate upon the square invertible 
case (which is much more complicated in these events) for k = 3 . 

Let B[, B 2 , ■. ■, B(, e M n (C) be given invertible matrices and consider the (non¬ 
linear) matrix equations 

1 . B\ — AiA 2 A3, 2 . B 2 = A 2 A3Ai, 3 . B 2 = A3AiA 2 , 

4 . B4 = AiA3A 2 , 5 . Z?5 = A 2 AiA3, 6. B(, = A3A 2 Aj 

in unknown matrices Ai, A 2 , A3 e M„( C). It is clearly necessary for a solution that 
Ai, A 2 and A3 be invertible. One natural generalization of Flanders’ question (which 
we consider) is the simultaneous solvability of all six equations. Of course, now 
with six data matrices and only three variable matrices, solvability should impose 
quite stringent conditions upon B \, B 2 ,..., B^. In particular, our assumption that 
the Bi' s (and therefore the A/’s) are square of the same size is necessary because of 
the requirements of conformability of the multiplication of the A/’s. We may also 
consider subsets of two or more of the six equations of the system ( 1 ). To this end, 
we first observe certain necessary conditions that flow from our equations. 

First of all, we have certain similarities among the B, ’s via A j ’s (which are as¬ 


sumed to exist): 

A { x B\A\=B 2 , A 2 l B 2 A 2 — B 2 , A 3 1 B 2 A 2 — Bi, 

Aj ‘b 4 Ai = B(j, A2 1 B$A 2 — B4, A^ 1 B(,A 2 — B5. 

This implies, in particular, that 

B[, B 2 and B 2 are mutually similar, 

i.e., have a common Jordan canonical form, ( 3 ) 

Z?4, R5 and B^ are mutually similar, 

i.e., have a common Jordan canonical form, ( 4 ) 

B\ and B4 are simultaneously similar to B 2 and B(, (via Ai), ( 5 ) 

B 2 and B5 are simultaneously similar to B 2 and B4 (via A 2 ), and (6) 

B 2 and B () are simultaneously similar to B\ and (via A3). ( 7 ) 


In fact, some 3 -tuple of matrices conveying the simultaneous similarities would have 
to give a simultaneous solution to Eqs. ( 1 ) or ( 2 ). 
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There are a number of other simultaneous similarity conditions involving various 
words in the B ,’s and their inverses, but all seem to be implied by (5)-(7). In view of 
(3) and (4), we may write 

B 2 = C 2 1 B 1 C 2 , B 3 = C 3 ' 5 iC 3 , B5 = C 5 1 B4C5, B(, = C 6 1 B^Cf,. 

( 8 ) 

Of course the Jordan canonical form of B\ and the one of B4 need not be the same 
(examples are easily given), but both Jordan canonical forms have the same determi¬ 
nant, as 

det Bj = det B /, all 1 ^ i, j ^ 6 . (9) 

All the necessary conditions mentioned thus far do not characterize the simulta¬ 
neous solvability of ( 1 ), as shown by the following. 


Example 1. Let 
B { = 


1/2 0 
0 2 


B 2 = 


0 

1/2 


B 3 = 


0 

1/2 


fl 4 = 


B 5 = 


b 6 = 


Then, det Bj = 1, all i , and 


( 10 ) 

(ID 


• B\ ~ B 2 and B 4 ~ B(, via 

• B 2 ~ Z ?3 and ZJ 5 ~ B 4 via 

• B 2 ~ B[ and B(, ~ B$ via 


• B 2 ~ Z ?3 and ZJ 5 ~ B 4 via the identity matrix and 

'0 r 
1 0 


In fact, up to scalar multiples, a calculation shows that these are the only similari¬ 
ties that work. Thus, all mentioned necessary conditions are met for a solution to (1). 


But the matrices 


, the identity matrix and 


are a commuting family, 


so that they (or any scalar multiples) cannot be a solution of ( 1 ). 

Moreover, if we consider the three matrices (10), we can see that the other three 
matrices ( 11 ) such that the system ( 1 ) is consistent, must have the form 


B 4 


a 0 
0 1 /a 


B5 


a 0 
0 1 /a 



0 

a 


However, these necessary conditions characterize simultaneous solvability for 
two natural types of subsets of our equations: transpositions (e.g., equations 1 and 4 
of (1)) and 3-cycles (e.g., equations 1-3 of (1)). The case of 3-cycles is essentially 
covered in [3], but we review it, with a complete description of all solutions, for 
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future use (observe that the simultaneous solution of the system ( 1 ) essentially in¬ 
volves coincidence of solutions for the two 3-cycles, 1, 2, 3 and 4, 5, 6 ). 

First of all, simultaneous solution of equations 1-3 of (1) is characterized by (3) 
and that of equations 4-6 of (1) by (4). For sufficiency, we note that 

Ai = BiC 2 , A 2 = C“ 1 C 3 and A 3 = C ” 1 ( 12 ) 

constitute a solution to equations 1-3 of (1), while 

A! = B 4 C 6 , A 2 = Cf 1 and A 3 = C~ l C 5 (13) 

constitute a solution to equations 4-6 of (1) (the matrices C 2 , C 3 , C 5 and C(, satisfy 
( 8 )). We note here that the situation is completely analogous for a single longer cycle. 

In fact, for the case of a 3-cycle, we may characterize all solutions. For an n-by-n 
matrix B, we denote by ( (i( B ) the group of invertible matrices that commute with B. 
Given two similar matrices A and B, denote the set of invertible matrices C such that 
C~ l AC = B by Z'(A, B), the set of matrices that convey the similarity. It is then a 
straightforward calculation that 

Lemma 2. If A, B e M n { C) are similar and Z e Z(A . B), then Z' e Z'(A, B) if 
and only if Z’ = XZY, in which X e @(A) and Y € ^?( B ). 

Proof. First, suppose that Z, Z’ e %{A, B). Then, from Z _1 AZ = B = (Z ') _1 AZ', 
we have Z’Z~ x A = AZ’Z~ { . Thus, Z'Z " 1 e #(A). We define X = Z'Z ” 1 and 
Y = I, to obtain Z' = XZY with X e ^(A) and Y e <#(B). 

Conversely, suppose that we can write Z' = XZY with X e r d(A) and Y e r d(B). 
Then, 

(Z'r 1 AZ' = (XZY)~ l A(XYZ) = Y~ 1 Z~ 1 X~ l AXZY 
= Y~ l Z~ l AZY = Y~ l BY = 5, 

as desired. □ 

We may then observe 

Lemma 3 

(a) If A\ = X\, A 2 = X 2 and A 3 = Xj constitute a simultaneous solution to equa¬ 
tions 1-3 o/(l), then Ai = Xi, A 2 = X 2 and A 3 = X 3 constitute another so¬ 
lution if and only if 

X x =T~ x XiT 2 , X 2 = Tf l X 2 T 3 and X 2 = T~ 1 X 3 T 1 , (14) 

in which Tj e ^( B , ), i — 1, 2, 3. 

(b) If A 1 = Fi, A 2 = Y 2 and A 3 = F 3 constitute a simultaneous solution to equa¬ 
tions 4-6 of { 1), then A\ = Y\, A 2 = Y 2 and A 3 = F 3 constitute another solu¬ 
tion if and only if 

F| = Tf'Yfl^, Y 2 = T~ l Y 2 T 4 and F 3 = T~%T 5 , 
in which T, e B ,), i — 4, 5, 6 . 


( 15 ) 
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Proof. We only prove (a); part (b) is similar. If X\, X 2 and X 3 form a solution to 
equations 1-3 and X\, X 2 and A 3 form another one, we can define 

T X = I, T 2 = X 2 X 3 X^Xp and T 3 = X 3 X~ l . 

It is easy to see that these matrices satisfy the equalities (14), and 7} e ( 6( Bf). 

On the other hand, a simple calculation shows that if the matrices Xj satisfy (14), 
Xt being a solution of ( 1 ), then they also are a solution of ( 1 ). □ 

The analog of Lemma 3 for «-cycles of square nonsingular matrices is also valid. 
Now, there is a solution of the system (1) if and only if there is an intersection 
among the triples that solve 1-3 and those that solve 4-6. If A; = Xj, i = 1,2, 3, 
is a particular solution to 1-3 and A,- = Yj, i = 1,2, 3 is a particular solution to 
4-6, then there is a common solution to 1-6 if and only if there exist T, e r S'(B ,), 
i = 1 ,..., 6 , such that 

Xi = Y u X 2 = Y 2 and X 3 = Y 3 , 

in which Xj, Yj, i = 1, 2, 3, are defined by (14) and (15). Letting .S', = Tf 1 , i = 
4, 5, 6 , algebraic manipulation shows that to solve the system (1) is equivalent to 
solving the system 

7iS 4 Fi =XiT 2 S 6 , T 2 S s Y 2 = X 2 T 3 S 4 , T 3 S 6 Y 3 = X 3 7)S 5 , (16) 

in which 7) e (li ,), i = 1, 2, 3, and .S', e (/l,), i = 4, 5, 6 , are unknown matrices. 
We return to this later to understand the solvability of system (1) in a particular 
case. 

We now know that a common Jordan form (i.e., mutual similarity) of the relevant 
Bj ’s characterizes the solvability of a cycle of equations (e.g., 1-3 or 4-6). This gen¬ 
eralizes in a simple way the square invertible case for k = 2 and may be generalized 
to the square invertible case for k > 3. What then about two of the equations, one 
from 1-3 and the other from 4-6. Each such pair amounts to transposition of the A 
indices. Somewhat surprisingly, condition (9) (for the relevant pair i, j) is sufficient 
(as well as necessary) for any particular pair of such equations. This may be seen 
using the observation of Sourour, which may be found in [5, p. 289]. 

Theorem 4. Suppose that A, B e M n (C) and that det A ^ 0. Then, for each non- 
cyclic permutation n I, there exist X \, X 2 , X 3 e M„(C) such that A = X \ X 2 X 3 
and B = X^i) 2Gr( 2 )2Lr(3) if and only if det B = det A. 


Proof. The necessity of the condition follows from the multiplicativity of the deter¬ 
minant, as mentioned in (9). For sufficiency, there are two cases to consider (a) one 
in which tt is the transposition 7 r( 3 ) = 1, tr(2) = 2 and tt( 1) = 3 and (b) the other 
in which n is an adjacent transposition. Since det A = det II f 0, det AT ? -1 = 1, 
we may write AB~ l = PQ [5, p. 289] in which P has the distinct nonzero eigen¬ 
values a\,a 2 , and Q has the eigenvalues l/oq, l/a 2 ,..., I/«„. Then Q~ l 

is similar to P, so that we may write S~ l PS = Q~ l for some S e M n (C). Letting 
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C = we have AC = P and SB = SQ~ l C~ x = PSC~' = AC(SC~ l ). 

Thus A = (SB)(CS~ i )C~ l and B = C~ l (CS~ l )(SB), which completes the proof 
in case (a). Case (b) may be verified by a modification of the same argument, but 
we mention another argument. Again, since det AS -1 = 1, we can write AB~ l = 
XY X~ x Y~ x for some X.Y e M n { C), is a multiplicative commutator if and only 
if its determinant is equal to 1 [5, p. 291]. Then Y~ l X~ l A = X~ l Y~ l B = some 
matrix Z, so that A = XYZ and B — YXZ, as was to be shown; the other adjacent 
transposition is similar. □ 

It is not difficult to see that there is an exactly analogous result for a longer string 
of X’s and a single transposition. Theorem 4 lies in contrast to Flanders’ observation 
for a transposition when there are only two matrices X\, X 2 (so that a transposition 
is also a cycle). Then the much stronger condition of similarity, in place of equal 
determinants, is required. 

We record as a theorem for comparison the case in which the permutation tc is a 
cycle. 

Theorem 5. The system of equations 1-3 has a solution if and only if B 1 , Ih and 
B 3 are mutually similar (i.e., have a common Jordan form). 

This theorem is an immediate consequence of Theorem 3.1 of [3], because the 
matrices If are invertible. 

Again, there is an analogous result for k equations corresponding to a k-cycle 
{B\ — A 1 A 2 ■ ■ ■ Ak, B 2 = A 2 -- AkA\, ..., Bk — AkA\ ■ ■ ■ Aa_i). For k — 3, The¬ 
orems 4 and 5 cover every possible pair of equations because every permutation is 
either a transposition or a cycle. For k ^ 4, this raises a natural question: what are 
the conditions on two (square, invertible) matrices A, B that they may be written 
A = X 1 X 2 ■ ■ ■ Xk and B — X- r (i ) X 7r ( 2 ) ■ ■ ■ X- T (k), in which it is a permutation that is 
neither a cycle nor a transposition. Always, det A = det It is necessary and A similar 
to B is sufficient. 

Corollary 6 . Suppose that A, B e M n (C) are invertible and n is a permutation of 
the set {1,2 ,,k}. Then A may be written A = X 1 X 2 ■ ■ ■ Xk while B is written as 
B — X 7T ^i)X 7T ( 2 ) ■ ■ ■ X 7I (k), in which X ,• e M,,( C), i = 1, 2,..., k, if and only if 

(1) B = A when n is the identity, or 

(2) B is similar to A when tt is a cyclic permutation (that is, 7 r(l) = i, 7 r( 2 ) = 
i + 1 , ..., jx{n) = i — 1 ), or 

(3) det B = det A, otherwise. 

We now turn to applying (16) to see more explicitly when the system (1) has a so¬ 
lution. Suppose that all of B\, ..., B(, are diagonal with distinct nonzero eigenvalues. 
Then we know that 
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B i = diag(A.i, ..., A„), 

B 2 = diag(A. CT( i),..., X a ( n) ), 

B 3 = diag(A. T( i), ...,X r(n) ), 

B 4 = diag(/xi, ..., ii n ), 

B 5 = diag(/ij(i), ..., fis ( n )). 

Bfr = diag(/XyQj,..., 

in which a , r, S and y are permutations and 

n n 

n =n as) 

;=i i=t 

(from equalities (8) and (9)). However, even in this special case, (8) and (9) do not 
insure a solution to (1). 

Let P a (respectively, l\, Pg, P Y ) be the permutation matrix corresponding to a 
(respectively, t, S, y); that is, P a = ]. Then, P a B\P'J = Ih ( P T B[ Pj = Bg, 

P S B 4 PJ = B S , PylhPy = B 6 ). 

Up to invertible diagonal multiples, these matrices provide the unique similarities 
in (8) because *^( B , ) is just the invertible diagonal matrices for each /?, . Now, with 
Xi = B\Pj, X 2 = Pn Pj- Xj = P z , Y\ = B 4 pJ, Y 2 = P S and T 3 = P y Pj (from 
(12) and (13)), we have, from (16), that Y\ and X\,Y 2 and X 2 , and F 3 and X 2 differ 
only by invertible diagonal multiplication, so that each pair has common support 
(zero-nonzero pattern). It follows that 

Pa = Py, Pa = P S Pr and P y = P T Pg, (19) 

in order for a solution to exist. 

In this case, we have the next important result. 

Lemma 7. Let B \. B 2 , Z? 3 , B 4 . B$ and B(, be defined as in (17), and assume that 
the conditions (19) hold. If (X\, X 2 , X 2 ) is a solution to equations 1-4 of the matrix 
system (1), then it is already a solution to the complete system. 

Proof. We know that (B\Pj, Pg, l\) is a solution to equations 1-3 of (1). So, ac¬ 
cording to Lemma 3, all solutions of these three equations have the form 

X, = Tf l B ] PjT 2 , X 2 = Tf ] PgTg and X 2 = T~ l P r T u (20) 

in which 7j, T 2 and 7 3 are nonsingular diagonal matrices. Clearly, 7) e ( Y> (If ) for 
all i = 1, 2, 3 and all 7 = 1,2,..., 6. 

Since (X 3 , X 2 , X 3 ) is also a solution to equation 4 of (1), we may write XiX 2 X 2 = 
B 4 . Then, 

X 2 X l X 2 = X 2 B 4 Xf 1 = t- 1 p s t 3 b 4 t- 1 pJt 2 
= T~ l P s B 4 PjT 2 = T~ l B 5 T 2 = B 5 
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and 


X 3 X 2 Xi = X; l B 4Xl = T-' P n B;'rB 4 T-' B ] P]T 1 
= T~ l P a B 4 P]T 2 = T~ l B 6 T 2 = B 6 , 

which completes the proof. □ 

Conditions (19) may be restated as 

Py = Pa, PaPs = P S Pa and P T = P n Pj, (21) 

a commutativity condition, together with a composition requirement. Thus, the sub¬ 
group of . ( /’ n generated by P a , P y , P$, P T is abelian and is generated by two of 
them, say P a , P$ (other pairs suffice). If this subgroup is decomposable, there will 
be further conditions on X\,..., X n and /i\. ..., /q, that refine (18) in order that 
(1) be solvable; if not, (17), (18) and (21) will be sufficient. In the former case, the 
additional conditions will be of the form 

]>=n^ 

in which each product is taken over indices i corresponding to an indecomposable 
component of the group generated by P a and P$. 

Proposition 8. Let B\, B 2 , If, If. B$ and B& be defined as in (17), and assume 
that the conditions (19) hold. If P a = Pa .l © Pa,2 and P$ — Pg.i © Ps.2, where P a j 
and Ps.i are nj-by-rij matrices (i = 1,2), then 

X[ ~'' X n ^ = [i\ ■ • • and = Bn\+\ *' * 

are necessary for the system (1) to be solvable. 

Proof. Suppose that (1) is solvable. Using (20), we have 

B 4 = (T ~ 1 B { PjT 2 )(T~ 1 P t T x )(T~ 1 P s T 3 ) 

= Tf ' If Pj 72 Tf 1 P a Pj 7, Tf 1 P, T 3 . (22) 

If we consider the partitions (remember that 7) and Bj are diagonal matrices) 

Th = Th.i © Th,i, h = 1,2,3 and Bj = Bj i © Bj j2 , j= 1.6, 

in which Th.i and B jj are «,-by-/z, matrices (i = 1, 2), the equality (22) splits into 
two equalities, corresponding to the blocks B 4 i and B 4 2 - From these equalities, it 
is straightforward to see that det B 4 j = det B\i,i = 1,2. □ 

In order to be explicit about the above remarks, we need to study when two per¬ 
mutations commute and, in particular, a joint decomposition involving cycles when 
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they do. We say that an n-by-n permutation matrix is irreducible if it corresponds to 
the permutation 

t 1 — ► *2 — ► ■ ■ ■ —► in —■► i t 

for i], (' 2 , ..., i„ some permutation of 1, 2,A cyclic permutation is a special 
case, and any irreducible permutation matrix is a permutation similarity of the stan¬ 
dard full cycle: i j = j, j = 1 , 2 ,..., n (which we typically assume it is, without 
loss of generality). First, a lemma is useful. 

Lemma 9. Let Q be a (0, 1 )-matrix which has at most one element equal to 1 in 
each row and in each column. If Q commutes with an irreducible permutation P, 
then either Q is a power of P or Q = O. 

Proof. It is well-known that any irreducible permutation matrix is nonderogatory. In 
this case, it is also known (see, for example, [4, p. 135]) that a matrix Q commutes 
with P if and only if the matrix Q is a polynomial in P of degree at most n — 1, 
where n is the order of P, 

Q = aol + a\P + • • • + a n —i P . (23) 

The structure of the matrices Q, I, P, P 2 , ..., P"~ 1 imposes that there is at most 
one nonzero coefficient a, in (23); moreover, the nonzero coefficient, if it exists, must 
be 1. □ 


Now, we may characterize a commuting pair of permutation matrices by showing 
that they may be put in a common block form in which commutation is obvious. 
The blocks of this form will correspond to the refined product conditions on the A.’s 
and pP s and, in turn, will lead to a characterization of the solvability of (1) for the 
problem in which the Bf s are as in (17). 

Theorem 10. Let P and Q be n-by-n permutation matrices. Then P and Q com¬ 
mute if and only if, up to simultaneous permutation similarity, 

P = P] ($...(£ P k an d Q = Q!®---®Q k , (24) 

in which P, and Q, are nj-by-tii permutation matrices {i = 1,..., k\ n = n;) 

and, for each i, one of the following possibilities holds'. 


(a) Pj is a standard full cycle and Qj is a power of P -,; 

(b) Pi = Ci ® • • ■ © C; and 

r oi i 


Qi = 



i 

o 


in which Cj is a standard full cycle and qi is a positive integer. 
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Proof. It is easy to see that if P and Q have the form (24) and each pair of matri¬ 
ces (P{ , Qj) satisfies either (a) or (b), then P and Q commute. We will prove the 
converse implication. 

Let T\ be a permutation matrix such that Pi = T\ P Tj = C ® P\ , where C is a 
standard full cycle. Let p be the order of C (we can identify C as having the greatest 
order among the standard full cycles in which the permutation represented by P 
decomposes, but it is not necessary). Let 

eff e ( i 2 r 

Qll Qll . ’ 

in which Q?) has the same size as C. 

The commutation condition, P\Q\ = Q\P\, implies that 

Cfiff = fiffc and P X Q^ = Q^P X . (25) 


Q i = t x qtJ = 


According to Lemma 9, we know that either Q } , = C q . for some q e Z + , or 

eff = o. 

If (j', V = C q , we take 


P 1 = C and Qi = C q . (26) 

In this case, it is clear that Q ( ^ and must be null matrices. Then, the matrices 

P| = Pi and Q\ — Qjj satisfy P\ = P\ ® P' v Qi = Qi ® Q[. Moreover, the size 
of the matrices P[ and Q\ is smaller than the size of P and Q and, according to (25), 
these new matrices are still commuting. 

Otherwise, if Q { ^ = O, it is possible to permute the columns of the matrix Q \ 
in order to obtain that the first p columns corresponding to the block Q\ t’ form the 
identity matrix. Observe that the permutation which is used in this reordering can 
be chosen so that it leaves the first p indices fixed. Thus, there exists a permutation 
matrix Tj such that 


P 2 = P 2 P,P' = 


and 


O 


Q 2 = T 2 QxV = 


Q 


o 



p(2) 

r 22 

p(2) 

r 2 3 


p( 2) 

r 32 

p( 2) 

^33 J 

0 

I 

0 " 

-,(2) 

t2\ 

0 

0 {2) 

2% 

0 

o (2) 
k?33 J 


The equality P2Q2 = Q2P2 implies that 


p( 2) 
r 22 


= c, 


CQ?) = Q\VC and P\t’Q 


■>( 2 ) 


(2) n (2) 


33 ^33 


n (2) p(2) 

^33 ^33 ’ 
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~ (2) 
Since P 2 is a permutation matrix, the first of these equalities provides that P 2 J = 

( 2 ) ( 2 ) ^ 

O and P 32 = O. Moreover, the second of them assures that either Q 2l — C q or 

( 2 ) n 

Q 2 y = O. In the first case, we get 


Pi = C © C and Q 1 


O / 

c q o 


If Q ( ~) = O, we repeat the process until we find a block ± O. It is clear 
that we will be able to find this nonnull block, since the successive matrices Q r are 
permutation (and finite) matrices. Thus, we will get 


O I 


Pi = C © • • • © C, Q 1 

-v-' 

r 



1 

o 


(27) 


and P and Q are simultaneously permutation similar to P\ © P[ and Q \ © Q\ , re¬ 
spectively, where the matrices P[ and O', commute and have smaller size than that 
of P and Q. 

Obviously, in every event, the matrices P and Q, up to simultaneous permutation 
similarity, can be written as 

P = P l ®P[, Q=Qi®Q\, 

in which P\, Qi are given either by (26) or (27) and the matrices P\ and Q\ com¬ 
mute. Then, we can repeat the above process on the pair of permutation matrices P\ 
and Q\ . Since each step of this process decreases the size of the considered matrices, 
the procedure is finite and completes the proof. □ 


Given a commuting pair of permutation matrices P and Q, according to Theorem 
10, there is a natural finest partition N\ U N2 U • • • U = N of N = {1, 2, ..., n} 
in which Nj is the set of indices corresponding to the blocks P, and Q, prior to the 
permutation of P and Q into the special form given in Theorem 10. We call this the 
commutation partition of P and Q, and, according to Proposition 8, this partition 
corresponds to the conditions on the k’s and /Ts needed to solve problem (1). Notice 
that if either P a or P$ is irreducible, the commutation partition is simply N, in which 
case the only condition will be (18), which is the successor to (9) in the case (17). To 
illustrate what can occur, we first give three examples. 

Example 11. Let 
"0 1 


1 0 


0 0 
1 0 
0 1 
0 0 


and P s = P; 
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Observe that the permutation P n is a standard full cycle (cr: 1—>2—>-3—>-4—* 1). 
According to (17) and (21), we have 


B i = diag(A.i, X 2 , A, 3 , X 4 ), 

B 2 = diag(A 2 , A 3 , X 4 , Xi), 

B 3 = diag(A. 4 , A-i, X 2 , X 3 ), 

B 4 = diag(/xi, M 2 , M3, IM), 

B 5 = diag(/X3, M4, Ml. M 2 ), 

B 6 = diag(M 2 , M3, M4, Ml)- 

We know that all solutions to equations 1-3 of (1) are given by (20). Let Tj = 
diag(f,i, ti 2 , f/ 3 , f; 4 ), for i = 1, 2, 3. 

The matrix equation X 1 X 3 X 2 = B 4 can be written as 

O 3 I 24 I 3 I _ Ml Il4l2ll32 _ M2 111 122^33 _ M3 *12*23*34 _ M4 

*11*23*34 M I 12 I 24 I 3 I ^ 2 ’ *13*21*32 A 3 ’ O 4 I 22 I 33 A 4 

(28) 

A simple calculation shows that the equality A. 1 A. 2 A 3 A .4 = M 1 M 2 M 3 M 4 i s the unique 
condition so that the system (28) be solvable. To obtain a solution to this system, we 
may take T\ — T 2 — / and 


T 3 = diag 


M4 M 1 M 2 M1M2M3 
ki X1A.2 Ai A2A.3 



Hence, the matrices 


"0 

0 

0 

M" 


■ 0 

0 

Ml M2 M 3 

A1A2A3 

0 " 

A2 

0 

0 

0 

, * 2 = 

0 

0 

0 

1 

0 

A.3 

0 

0 

Ml 

h 

0 

0 

0 

0 

0 

X4 

0 


_ 0 

MlM 2 

L1X2 

0 

0 


*3 


0 

M1M2 

0 

0 


0 0 

0 0 

AlAgAl Q 

m 1/22/23 

0 1 


Ai" 

Ml 

0 

0 

0 


form a solution to equations I -4 of (1). Lemma 7 assures that this solution is already 
a solution to the complete matrix system (1). 
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Example 12 . Let 


"0 

1 

0 

o' 


'1 

0 

0 

0" 

1 

0 

0 

0 

and Pg = 

0 

1 

0 

0 

0 

0 

0 

1 


0 

0 

0 

1 

0 

0 

1 

0 


0 

0 

1 

0 


Observe that the permutation matrices P a and Pg are decomposable: P a = C © C 
and Pg = C 2 © C, where C is the standard full cycle of order 2. In this case, we have 

B[ = diag(Lj, X 2 , A. 3 , X 4 ), 

B 2 = diag(>.2, A,i, A4, X 2 ), 

Bj = diag(L 2 , A.i, X 3 , X 4 ), 

B 4 = diag(/xi, m 2 , /x 3 , n 4 ), 

B 5 = diag(/xi, H 2 , M4, M 3 ), 

= diag(M2, mi, M4, M3)- 


Again, all solutions to equations 1-3 are given by (20). Now, the matrix equation 
X 1 X 3 X 2 = B 4 implies that 


122*31 _ Ml *21*32 _ M2 *14*33 _ M3 *13*34 _ M4 

*21*32 M *22*31 ^2 *13*34 ^3 *14*33 ^4 

Notice that to solve this system, the conditions X\X 2 = M 1 M 2 and A 3 L 4 = M3M4 are 
necessary. The matrices T\ — T 2 — I and 


t /Ml . M3 . 

73 = d,a Hv *’ I? 1 

provide a solution to the matrix system: 


"0 

^1 

0 

0" 


r/M 

Ai 

0 

0 

0" 

X 2 

0 

0 

0 

, x 2 = 

0 

1 

0 

0 

0 

0 

0 

^3 

0 

0 

0 

1 

0 

0 

X 4 

0 


0 

0 

m 

^3 

0 


0^-00 

Ml 

10 0 0 

0 0^0 

M3 

0 0 0 1 


Example 13 . Finally, let 


"0 

1 

0 

o ' 


'0 

0 

1 

0 " 

1 

0 

0 

0 

and Pg = 

0 

0 

0 

1 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 


0 

1 

0 

0 
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(these matrices have the structure given in (27), where C is the standard full cycle of 
order 2 and q — 2). Then, 


B i = diag(Ai, A. 2 , 7-3, ^ 4 ), 
#2 = diag( 7 . 2 , A.i, A 4 , 13 ), 
Z ?3 = diag(l 4 , A. 3 , X 2 , M), 

B 4 = diag(/ri, M 2 , M3, M4), 
B 5 = diag(/43, M4, Ml, M 2 ), 
= dhig(M2, mi, M4, M3>- 


The matrix equation X 1 X 3 X 2 = B 4 can be written as 

03^22^31 _ Ml M4?21?32 _ M2 01^33 _ M3 712*23*34 _ M4 

*11*23*32 ^-1 *12*24*31 ^-2 *13*21*34 7-3 *14*22*33 ^4 

This system is solvable if and only if A 1 A 2 A 37.4 = M 1 M 2 M 3 M 4 - A solution to it is 
given by 


7| = diag 



AlT-2 

M 1 M 2 



and 


T 2 = Id 


r 3 = diag 



I1A.2I4 

M1M2M4 


Then, the matrices 


*1 


0 Ai 0 0 

® 0 0 0 

0 0 0 a .3 

0 0 A 4 0 


*3 


"00 0 r 
00^-0 
0^00 
10 0 0 


*2 


"0 0 
0 0 
1 0 


*4*,2^4 q- 

^ 1 / 42/44 

0 1 

0 0 

0 0 


form a solution to the matrix system ( 1 ). 


Our aim is to find a necessary and suficient condition so that the system (1) be 
solvable when the matrices B , are as in (17). We know that conditions (19) are nec¬ 
essary. From these conditions, the permutation matrices P a and Pg commute. Then, 
according to Theorem 10, there exists a permutation similarity S such that 

S~ l P a S = P a — P a , 1 ® ■ ■ • © P a ,k and S~ l PgS = P s = P s j © • ■ ■ © Pg^, 
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where P a j and Pgj are iij-by-n, permutation matrices, for i = \..... k, and each 
pair ( P a i , PsP has either the form (C, C q ) or the form given in (27); in both cases, 
C is a standard full cycle and q is a positive integer. 

If we define Bj — S ~ 1 If S. for i = 1, ..., 6 , and (X \. X 2 , X 3 ) is a solution to the 
system 

B l =X l X 2 X 3 , Bi = X 2X3X1, B 3 = X 3 X l X 2 , 

B a = X x X 3 X 2 , B 5 = X 2 XiX 3 , B 6 = X 3 X 2 X l , 

then X[ = SX[S~ l , X 2 = SX 2 S~ l and X 3 = SX 3 S' -1 form a solution to (1). More¬ 
over, notice that B 2 — P a B[Pj if and only if B 2 = P a B[Pj, and so on for the rest of 
the known relations involving the Bj ’s. Therefore, without loss of generality, we may 
suppose that P a = P a and Pg = Pg. Then the system (1) splits into the subsystems 

B U = X u X 2 J X Xi , B 2J = X 2 J X 3 J X Ui , B Xi = X 3 J X u X 2 ,i, 

B 4 ,i = X u X Xi X 2<i , B s ,i = X 2 <i X hi X Xi , B 6J = X 3 jX 2 jX hi 

for i = 1-- k, in which B 2 j = P a y B^ Pj t , B 3 j = P r jB u Pjj, B 5 j = 

Ps.iB 4 .iPjj and « 6 , ; = P 0 .iB 4 .iPjj. 

So, we only need to study the solvability of (1) in two cases: 

(1) P„ is a standard full cycle and Pg is a power of P a , and 

(2) the pair ( P a , Pg) has the structure given in (27). 


Notice that (22) can also be written, after some algebraic manipulations, as 
B 4 B;' = Tf\pjT x Pg){PjT 2 P a ){PjT~ x Ps)(PjT~ l P a )T 3 . 

From Examples 11-13, it appears possible to consider 73 = /; we will test this 
idea. Under this supposition, the equation to solve is 

B 4 BJ 1 = T-\pjT { Pg)(PjT~ { P a )T 3 , 
which is equivalent to 


1 (i) ' U,i 
h,i ■ h.o-Hi) 


— for i = 1 ,... ,n. 


(29) 


Proposition 14. Let B 1 , B 2 , B 3 , B 4 , Bg and Bg be defined as in (17) and assume 
that the conditions (19) hold. Assume also that P a corresponds to a standard full 
cycle of order n and Pg is a power of P a . Then , the system (1) is solvable if and only 
if det B 1 = det By. moreover, a solution to the system is given by 


X 1 — B\Pj, X 2 = P S T 3 and 


where 


r3 = diag (l? 


Ml M2 
klT-2 


Ml ■■ 'Mn-l 
A.t • • • A. n _i 


II 

oTj 

(30) 

7' O' 

(31) 
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Proof. The condition det B\ — det B 4 is obviously necessary. Moreover, using the 
matrices defined by (30), it is straightforward to see that 

X l X 2 X 3 = B l , X 2 X 3 X l = B 2 and X 3 X 3 X 2 = B 3 , 


and 


* 1 * 3*2 = B.P^T-'P t P & T 3 = B { PjT~ x P a T 3 . 


(32) 


Then, the definition of T 3 given in (31) and the special form of P a (that is, er(i') = 
i + 1 for i — 1 ,..., n — 1 and a in) = 1 ) provide 


PjT~ l P a = diag ( 1, —,. 

Ml 


*i 


*1 • • • *71-2 *1 • • • *71-1 


Ml ‘ ‘ * Mil — 2 Ml * * ‘ Mw— 1 


(33) 


From (17) and (31)-(33), 

*1*3*2 = diag ^Mi> M2, • • ■, Mn-1> 
= diag ( Ml, M2, • ■ •, Mn-i- 


*1 *17 

Ml •• • Mn-l 
M ,7 det B\ 
det B 4 


Then, det B\ = det If is a sufficient condition so that X\ X 3 X 2 = B 4 . Lemma 7 
completes the proof. □ 


To find a solution to system (1) in the case in which P a and P$ have the structure 
given in (27), we need some preliminary results. 

As usual, we denote by |_vj the integer part of x and by mod(x, y) the remainder 
of x after division by y. 


Lemma 15. Let P a = C ® • • • ® C , in which C is a standard full cycle of order p 

" v “ ^ 

r 

and let 1^/5) pr. Then , 


(!) i = L '-y\P + mod (' - 1, P) + 1> 

(2) er -1 (0 = \rj-\p + mod); - 2, p) + 1. 


Proof. Part (1) is straightforward to prove: 

i — 1 


i — (i — 1 ) + 1 — 


P J 


p + mod(i — 1 , p) + 1 . 


To prove part (2), observe that the permutation a may be described as follows: 


1 


2 p —► 1, (p + 1) -* (p + 2) • • ■ -* 2p (p + 1), 

[(r - 1 ) p + 1 ] -> [(r - \)p + 2 ] -* ■ • ■ -* rp -> [(r - 1 )p + 1 ]. 


We analyze three cases. 
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• If i = ap + 1 with 0 ^ a < r, we have a 1 (i) = i + p — 1 and 

i — 1 

- p + mod(/ — 2, p) + 1 = ap + (p — 1) + 1 = i + p — 1. 

_ P J 

• If i — ap + h with 0 ^ a < r and 2 ^ h ^ p — 1, we know that er -1 (i) = i — 1; 
on the other hand, 


i — 1 


p + mod(; — 2, p) + 1 = ap + (h — 2) + 1 = ap + h — 1 = i — 1. 


P J 

• If i — ap with 1 ^ a ^ r, then a -1 (i) = i — 1 and 
i — 1 

p + mod(; - 2, p) + 1 = (a - 1 )p + (p - 2) + 1 
= ap — 1 = i — 1. □ 


L P 


Lemma 16. If p,r, i e Z + , we can write 
i = (y - 1 )p + modO -k,p) + 1, 

in which 


(34) 


J 


i — 1 
L P J 


+ 1 and k = mod ^ 1 + 


i — 1 


- i, P + 1. 


Proof. According to part (1) of Lemma 15, if we write i as proposed in (34), we 
have 


7-1 = 


i — 1 


and modO — k, p) = mod(i — 1, p). 


The second of these equalities implies that mod(£ — 1, p) = mod( / — i, p), that is, 

i — 1 


mod(& — 1, p) — mod I 1 + 


P J 


~t,P 


In order that this relation be true, we may take 


k = mod(l + L^y-J - i, p) + 1. 


□ 


Fixing p e Z + , we define the functions J p (i) and K p {i) as 


Jp(i) 


! — 1 


P J 

They are the unique values such that 


+ 1, Kp(i) — mod ( 1 + 


i — 1 
L P J 


-i, P 


i = [J p (i) - 1 ]p + mod[7 ; ,(0 - K p (i), p]+ 1, 
JpU) > 1 and U K p (i) < p. 


( 35 ) 


It is easy to see that if i = (j — 1 )p + mod (j — k, p) + 1, then 
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a J (i) = O' - 1 )p + modO - k - 1, p) + l. 

In other words, 

Jp(a~ l (0) = Jpti) and K p (a~ l (i)) = mod(K p (i), p) + 1. 

Now, we need to study the permutation 8 . Let 8 be the permutation corresponding 
to C q . Then, <5 can be described by 

1 -►(/>+ 1) -► (2p + 1) [(r - l)p + 1] -> 5(1), 

2 -* (p + 2 ) (2p + 2 ) -»-> [(r - Dp + 2] -> 5(2), 

p -+ 2p -»• 3p - > rp -* 8 (p), 

where, obviously, 1 ^ S(i) f p, tor i = 1,2,..., p. From these chains of images we 
may construct the cycles of <5, according to the values 8 (i). It is easy to see that if <) is 
an irreducible permutation, then S is also irreducible; in this case, the solution of the 
matrix system studied can be constructed as in the case solved above, interchanging 
there the roles of T\ and T 3 . 

For p < i ^ rp, it is straightforward to prove that J p ( 8 ~ l (D) = J P (D — 1 and 
K p ( 8 ~ l (i)) = mod (K p (i) - 2, p) + 1. 

To simplify the notation, we define the values, for n e Z + and m = 1,2,...,/?, 

, , k.(n-l)p+mod(n-m,p)+l 

a(n,m) = -. (36) 

l^(n — l)/?+mod(«— m,p)-\-l 


Lemma 17. The values defined in (36) satisfy the following properties : 

(1) K p (D ) = it; 

( 2 ) n:=i ot(n, m) = U n „L(n-i )p +i 

Proof. Part (1) follows from (35). To prove part (2), notice that, for m — 1,2,...,/?, 
(n — Dp + 1 ^ (n — Dp + mod(?? —m,p) + 1 ^ np. 

Since mod(n — m\, p) ^ mod(?7 — w 2 , p) when m\ ^ m 2 , the expression given by 
(n — Dp + mod(n — m, p) + 1 

takes p distinct values, that is, all the integer values from (n — 1)/? + 1 to np inclu¬ 
sive. □ 

Proposition 18. Let Bi, IL. B 3 , B 4 , B 5 and B(, be defined as in (17) and assume 
that the conditions (19) hold. Assume also that P a and P$ have the structure de¬ 
scribed in (27), in which C is a standard full cycle of order p. Then , the system (1) 
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is solvable if and only if det B\ — det B^\ moreover, a solution to the system is given 
by 

Xi = Tf'BxPj, X 2 = P S T 3 and X 3 = Tf'P t T u 


where Th — diag(?/,,i, th,2, ■ ■ ■, tk,rp), h = 1, 3, and 


t\,i — 


h,i = 


■Ip(0-P x 

— if K p {i) = p and 1 < J p (i) 

s= 1 /A 

1 otherwise, 

1 


\Jp(i)-\\p 

n 


5=1 



fK p (i )-1 


a(J p (i), m) 


m= 1 


<r-l. 


if K p (i) = 1, 


if 2 < ^p(i) < p. 


(37) 


(38) 


Proof. To prove this, we determine a sufficient condition so that the values given by 
(37) and (38) satisfy Eqs. (29). 


(1) For 1 ^ i ^ p, we know that J p {i) = 1 and (r — l)p + 1 ^ <5 1 (;') ^ rp. Since 




S~Hi)~ 1 

P 


+ 1 = r. 


in the definition (37) we find that fj = 1. We need still to study three sub¬ 
cases. 

(a) If K p (i) = 1, we know that i = 1. Since K p (c r —1 (1)) = 2, according to 
(38), we have 

^1,5-!(1) ' r 3,1 _ 11 _ Ml 

f U ' f 3,o-!(i) 1 • a{J p {\), 1) A.i 

(b) For 2 ^ K p (i) ^ p — 1 (that is, for 3 ^ i ^ p), we know that 
K p (a~ l (i)) = Kpii) + 1. Then, 

^l.g-TO • ? 3,1 _ 1 ■ nf=i° ^(TpOQ.tw) _ 1 _ 

h,i ■ t 3 ,o-fi) ~ 1 ■ n K J? 0. rn) _ ^/> (/)) _ 

(c) If Kp(i) = p, that is, if i = 2, we have K p (ct~ l (/)) = 1. Then, 

f l,3~ 1 (2) • r 3.2 _ 1 • n»i = l a ( J p( 2 )’ m ) _ Ylm=\ a ( J p( 2 )’ m ) 

■ ^,«r-i(2) “ nf=i t • i ” «(4( 2 ). p) ■ nf=i k 

- 1 _ M2 

a(J p (2), p) \ 2 
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(2) For p + 1 ^ i ^ rp, we know that 2 ^ J p (i ) ^ r and J p (8~ l (i )) = J p (i ) — 1. 
Then, 

(a) If K p (i) = 1, we have K p (er _1 (/)) = 2 and K p (8~ l (i)) = p. Thus, 


r i,6- 1 (0 ' f 3,i 


i-r[^p(0 1]P Aj_ , 

1 ii=l /-t.v ' 1 


(rti“ ,_11 " t) “wo. i) 


fM 

V 


(b) For 2 ^ K p (i) ^ p — 1, we know that /^(cr '(i)) = K p {i) + 1 and 
KpiS-Hi)) = Kp(i)~ l.Then, 


f l,6- 1 (0 ' f 3,/ 

f i,i ‘ 


(r£P~'«wo.»'>) 


o(/ P (i), ^(0) 


/ M_ 
Xi' 


(c) Finally, we study the subcase in which K p (i) — p. Since K p (cr 1 (i)) = 1 
and K p ( 8 ~ 1 (i )) = /? — 1, we can write 


r i,6- 1 (0 ‘ f 3,i 
f U ’ ^ 3 ,<t 1 (i) 


(n 


[/ P (i)-i]p i, 

J=1 u. 


)(n 


/! = i a(J p (i), m)) 


t\,i * 1 


r-r^p(') A*_ 

1 jJL s 

t\,i ■ <x(J p (i), p) 


In this expression, we need yet to consider two cases: J p (i) < r and J p (i) = 
r. When J p (i) < r, we have 

ks_ r-f Jpd) Aj_ 

1 l,y=l jji s _ 1 l.v=l p, s _ 1 _ fJi 

t Ui -a(J p (i), p) (n/li 0 a(J p (i), P) <x(Jp(i)’K p (i)) X,' 

On the other hand, when J p (i ) = r, we know that t\j = 1. Then, 

t i,8~Hi)' t 3’ i _ _ detfil ^ 

l U ' h.u-HO a(J p (i), K p (i)) det Z? 4 X t ' 

This is the unique case in which we need an additional condition so that the 
corresponding equation given in (29) be satisfied: it must be that det B\ = 
det Z? 4 . 


Hence, it is clear that the necessary condition det B\ = det B 4 is also a sufficient 
condition for solvability of (1). □ 
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Using the results found above, we have proved our final main theorem. 

Theorem 19. Let If. ..., If , be as in (17), with '/. \ , ... , k n distinct and /i] ,..., /i n 
distinct, and let P a , If, P$, P y be the permutation matrices defined by Bi, B$, If 
and If, (i. e ., If — P a B\Pj, If = P T B\Pj , If = Pi, If Pj and If, = P y If Py). 
Then , the problem (1) has a solution if and only if 

(1) P a = P Y — If Pi, = If If and 

(2) Q j eN . kj = ri/siV, P j for each N, in the commutation partition of the commut¬ 
ing pair of matrices P (j , P/,- 

Finally, we note that there is a simple corresponding result if B\,...,B(, 
are simultaneously diagonalizable by similarity and the diagonal forms meet the 
conditions of Theorem 19. 
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